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1.0  SUMMARY 


The  Two-Dimensional  Douglas  Neumann  Program  for  calculating  the  potential 
flow  about  bodies  of  arbitrary  shape  has  been  extended  to  handle  lifting,  in¬ 
finite  cascades.  The  resulting  very  general  program  allows  a  wide  range  of 
heretofore  intractable  problems  to  be  solved.  Essentially,  the  program  can 
handle  any  problem  in  which  the  flow  pattern  repeats  indefinitely  along  an  axis. 
In  the  Douglas  Neumann  program  this  axis  is  the  y-axis.  This  permits  the  calcu¬ 
lation  of  the  flow  about  a  lifting  or  nonlifting  cascade  having  any  stagger  angle 
and  spaci  lg  and  having  arbitrary  blade  geometry.  The  program  can  also  calculate 
the  flow  about  more  than  one  cascade.  Thus  it  can  handle  ihe  interaction  prob¬ 
lem  of  two  or  more  parallel  cascades. 

The  Douglas  Cascade  Program  is  compared  with  other  theoretical  methods, 
special  analytical  cases,  and  experimental  data.  Program  details,  among  which 
are  input-output  format  and  FORTRAN  listing  are  given  in  the  appendices. 
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4.0  NOTATION 


'-o 

C1 

c2 


i,  J 

i 

K 

M 

n,  t 
NC 

q. 

s 

SP 

U 

UI 

u„ 


Kutta-condition  matrix  defined  in  appendix  A 

chord  length  of  cascade  blade 

complex  coordinate  of  body  element  midpoint 

complex  coordinate  of  body  element  first  endpoint 

complex  coordinate  of  body  element  second  endpoint 

lift  coefficient  per  cascade  blade  normalized  with  U  the  average 

2_r 

velocity  modulus,  yc  .  If  some  other  normalizing  velocity, 

aru  c 

is  used  the  lift  coefficient  is  ,  .  When  a  set  of  cascades  is 

TJ*C 

considered,  is  the  lift  coefficient  for  the  set. 

pressure  coefficient  having  the  average  onset-flow  velocity  as  the 

normalizing  velocity 

unit  vectors  in  the  x-and  y-directions,  respectively 
complex  unit, 

midpoint  and  element  subscripts,  respectively 
complex  source  strength,  K  =  M  +  iT 
source  strength  for  point  source 

unit  vectors  in  the  normal  and  tangential  directions,  respectively 
number  of  cascades  considered 
surface  location  of  a  source 
surface  location  of  a  general  point 

cascade  spacing,  the  flow  pattern  repeats  with  this  spacing  along 
a  prescribed  axis 
modulus  of  V 
modulus  of  Vj 

modulus  of  any  normalizing  velocity  V_ 
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V(q,s) 


z 

a 

“e 

°k 

°r 

Ax 

V 

AVO 
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AV90 

IB 


AY 

1 

* 
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X 


average  onset-flow  velocity 


Y  +  V 
I  E 


velocity  at  s  due  to  a  source  at  q 

exit  velocity,  velocity  at  x  ■  go 

inlet  velocity,  velocity  at  x  *  -co 

invluenee  of  element  k  at  point  a 

influence  of  element  k  at  the  Midpoint  of  eleaent 


i, 


coordinates  of  a  general  point 
x  +  iy,  complex  coordinate  of  a  general  point 

average  angle  of  attack  measured  in  +v*e  counterclockwise  direction 

exit  angle  of  attack 

angle  of  attack  of  the  kth  eleaent 

inlet  angle  of  attack 

cascade  turning  angle 

circulation  about  an  individual  cascade  blade 

effect  of  a  uniform  onset  flow  at  zero  angle  of  attack  on  the 

Kutta  condition  of  the  mth  cascade  in  a  system  of  several 

cascades 

o 

effect  of  a  uniform  onset  flow  at  90  angle  of  attack  on  the 
Kutta  condition  of  the  mth  cascade  in  a  system  of  several  cascades 
effect  on  the  Kutta  condition  of  the  mth  cascade  due  to  circula¬ 
tory  flow  about  the  nth  cascade  in  a  system  of  several  cascades 
£  +  ir£  ,  complex  coordinate  of  a  source 

stagger  angle  (see  figure  3)>  measured  clockwise 
surface  source  distribution 
airfoil  trailing-edge  angle 
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5.0  INTRODUCTION 


The  Douglas  Neumann  Program  is  a  powerful  tool  for  determining  the  potential 
flow  about  one  or  more  arbitrary  two-dimensional  lifting  bodies.  The  Neumann 
program  is  rigorous  in  the  sense  that  the  exact  solution  is  approached  in  the 
limit  as  the  number  of  points  describing  the  body  goes  to  infinity.  This  power¬ 
ful  method  has  now  been  extended  to  calculate  the  flow  about  infinite  two-dimen¬ 
sional  lifting  cascades.  Cascade  data  are  ued  to  advantage  when  working  with 
compressor  stages,  turning  vanes,  or  propeller  blades.  A.  J.  Acosta  (Calif. 

Inst,  of  Tech.)  and  H.  P.  Linhardt  (reference  l)  have  concluded  that  the  use 
of  cascade  theory  to  predict  propeller  characteristics  is  accurate  even  when 
the  propeller  is  extreme  in  configuration.  These  authors  tested  a  low-aspect- 
ratio  axial-flow  pump  propeller.  Comparison  of  the  experimental  results  with 
three-dimensional  vortex  and  two-dimensional  cascade  theory  showed  the  simple 
cascade  theory  to  be  superior. 

A  cascade  is  defined  as  simply  a  series  of  identical  bodies  equally  spaced 
and  identically  oriented.  There  are  no  restrictions  on  the  body,  shape,  spacing, 
and  stagger  angle.  The  cascade  of  bodies  may  be  lifting  or  nonlifting.  In 
general,  the  program  can  handle  any  problem  in  which  the  flow  pattern  repeats 
along  an  axis  from  plus  to  minus  infinity.  Figure  3  shows  a  double  cascade 
along  with  a  graphical  representation  of  the  cascade  parameters. 

Shown  in  figures  4,  5,  6,  and  7  are  some  of  the  extreme  configurations  a 
cascade  might  take,  given  the  flexibility  this  program  affords.  The  program 
can  also  be  used  to  calculate  the  flow  past  a  series  of  cascade  stages,  i.e., 
more  than  one  cascade.  The  stages,  however,  cannot  move  relative  to  each 
other,  since  this  would  be  an  unsteady- flow  problem. 

In  some  cases  the  effect  of  boundary- layer  displacement  thickness  on  the 

cascade  blade  is  important.  In  these  cases  the  displacement  thickness  can  be 
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added  to  the  blade  and  the  resulting  thicker  blade  can  be  used  in  the  cascade 
program.  This  displacement-thickness  technique  has  been  tried  successfully 
on  a  single  isolated  airfoil. 

Hie  Douglas  levaann  program  vith  its  cascade  modification  calculates  the 
velocity  and  pressure  distribution  normalized  to  average  velocity  and  lift 
coefficient  and  moment  coefficient  per  cascade  blade.  Also  calculated  are  the 
inlet  and  exit  velocities  and  the  cascade  turning  angle.  Appendix  C  gives  the 
FQRTRAlf  listing  and  all  the  details  of  the  input  and  output  of  the  Douglas 
Cascade  Program. 
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6.0  THEORY 


The  technique  employed  by  the  Neumann  Program  to  solve  the  fluid-flow 
problem  is  to  apply  a  source  distribution  of  appropriate  strength  cn  the  sur¬ 
face  of  the  body  in  such  a  vay  that  the  flow  normal  to  the  surface  of  the  body 
is  either  zero  or  prescribed.  This  technique  is  described  in  great  detail  in 
references  2,  3,  and  4.  When  the  Neumann  boundary  condition  is  applied,  an 
integral  equation  in  source  strength  C"  is  obtained.  This  integral  equation 
is 


Vo  •  n  -  C5-(S)  +  /erf?)  A(?,s)  d 


cn 


BODV 


where  A(q,s)  =*  n*  V(q,s)  and  Vqq  is  the  onset  flow. 


In  the  unmodified  program  V(a,s)  is  the  familiar  velocity  at  s  due  to 
a  unit  source  at  q.  If  x,y  are  coordinates  associated  with  s  and  ^  ? 
with  q,  then 

V(g,s)=  tVx+IvY 


where 


v*  = 


vY  = 


(*■-$) 


(2) 


(x-5)2+(y-tjy 


or,  in  complex  notation. 


V  -  L  V  - 
'x  L  VY  ~ 


To  modify  the  program  to  handle  infinite  cascades  a  new  velocity  at  a  point  s 


due  to  a  unit  source  at  q  in  a  cascade  is  used. 


This  "new"  velocity  is  the  sum  of  the  velocities  due  to  a  row  of  sources 
equally  spaced  along  the  y  axis.  The  sum  is  tne  series  representation  of 
the  hyperbolic  cotangent  function  (see  Lamb,  reference  ^ ,  page  7l).  Thus  the 
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cascade  velocity  is 


V*  '^Y  “  2S p  i  SP  Z(s)-%  (tj)  f 


1 

2SP 

3  (*• 

■  %)]  aLyiM.  [ 

.SP  V*  1 

1 

M 

i _ i 

H>] 

F* 

j-<UAb/L 

[&(*-§)] 

F 

I 

2.SP 

■*rf  , 

_sp  (!/ 

:)’■ 

« 

AJUYL 

JL  , 

SP  1 

(y-D. 

< 

r 

Aim, 

>(n) 

2 

-  -f 

SP 

i 

Afr 

Is  i 

-  2 

► 

J 

where  SP  is  the  cascade  spacing.  The  cascade  vortex  velocity  is  (Lamb,  page 
244)  simply  the  cascade  source  velocity  rotated  90°  clockwise. 


The  technique  employed  to  solve  the  integral  equation  given  in  (l)  is  to 
approximate  the  body  surface  by  straight-line  elements.  The  source  strength  is 
assumed  constant  along  any  one  element,  but  it  varies  from  element  to  element. 
Now  if  the  Neumann  boundary  condition  is  applied  at  the  midpoints  of  each  of 
these  elements,  equation  (l)  can  be  written  as 


Here  we  note  that  A(q, s)  is  now  written  Aj(q)  since  the  positions  of 
the  general  points  are  now  fixed  as  the  element  mid-points,,  Referring  to  the 
definition  of  A(q,s)  given  for  (l),  we  write 


Aj  l<fc)=  n  •  Vj  (a.) 


(5) 


V.(q)  is  the  velocity  at  the  mid-point  of  the  jth  element  due  to  a  unit  source 
J 


at  q , 
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Let  the  quantity  Wjk  be  defined  as 


V  7 


lx  +  Ly 


(6) 


For  convenience  the  complex  form  of  W(q,s  )  or  W j ( q )  and  V(q,s  )  or  Vj(q) 
will  be  adopted  and  used  henceforth.  Substituting  the  cascade- velocity  source 
function  found  in  (3)  into  (6),  we  have  the  following 


The  kth  element  over  wr.icr.  tre  integration  is  to  be  performed  is  shown  in 
figure  1.  F.-om  figure  1  the  following  relations  are  evident 


Figure  1.  -  A  typical  straight  line  element  of  the  body  surface. 
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If  the  relations  of  (8)  are  used,  equation  (?)  may  be  rewritten  as  follows: 


W  =  e  k  j  -r_  qjMi 


’2k 


25P 


2L 

SP 


which  upon  integration  becomes 


w  =  X-  ZY  = 


-toe, 


2^ 


Yn  < 


~ 

_SP 

(zj  ~c,n)" 

LI 

AUYlJb 

n—  1 
_  SP  ' 

:zrcuJ 

II 

(9) 


Equation  (9)  is  the  basic  cascade  source  function  used  in  the  Douglas  Cascade 
Program.  It  represents  the  complex  velocity  at  the  jth  element  midpoint  Zj, 
due  to  the  kth  source  element  in  cascade. 


If  we  now  use  the  definition  'iven  in  equation  (6),  equation  (4)  can  be 
rewritten  thus: 


N 


k  =  1 


Ni 


(10) 


Equation  (10)  is  then  solved  using  equation  (9)  for  the  unicnown  cr^.  Once 
theCT.  values  are  known,  tne  velocity  and  pressure  anywhere  in  the  flow  field 
can  be  calculated. 


To  solve  the  general  case  of  a  lifting  cascade  at  any  angle  of  attack, 
"basic"  flows  are  calculated  and  superimposed  in  such  a  way  that  the  correct 
angle  of  attack  is  obtained  and  the  Kutta  condition  is  satisfied.  These  "basic" 
flows  are  the  following: 

1)  Flow  at  zero  angle  of  attack. 

2)  Flow  at  90°  angle  of  attack. 

3)  Circulatory  flow  for  each  cascade. 
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Equation  (10)  is  solved  for  each  one  of  the  basic  onset  flows;  here  is  the 
onset  flow  in  the  equation.  Superposition  of  solutions  is  possible  because 
the  potential  equation  is  linear  and  the  boundary  condition  on  the  cascade 
blades  is  homogeneous.  The  details  of  the  superposition  technique  are  given 
in  Appendix  A. 


In  a  cascade  there  is  an  infinite  number  of  airfoils  or  blades,  each  hav¬ 
ing  a  circulation.  Therefore  since  the  cascade  runs  along  the  y-axis,  there 
exists  an  upwash  an  infinite  distance  upstream  of  the  lifting  cascade  and  a 
dovnvash  an  infinite  distance  downstream.  At  these  distances,  tne  lifting 
airfoils  act  like  a  row  of  equally  spaced  vortices.  The  magnitude  of  the  upwash 
and  downwash  due  to  this  row  of  vortices  can  be  deduced  from  Lamb  (reference  5) 
as 


^up  —  ^cLotirn.  ~ 


2SP 


Uri  Cl^c 
4-U  SP 


(11) 


where  P  is  the  circulation  per  cascade  body  and  SP  is  the  cascade  spacing. 

Uy^is  any  convenient  normalizing  velocity.  In  the  Douglas  program  U  is  Just 
L  NC 

u  .  (When  more  than  one  cascade  is  involved,  T  is  replaced  by  ^  pm;  KC 


m=l 

is  the  number  of  cascades.  Also  is  the  lift  coefficient  for  the  set  of 
cascades.)  The  cascade  therefore  turns  the  flow.  The  upwash  and  downwash  and 
the  cascade  turning  angle  are  shown  in  figure  3.  For  a  cascade,  the  lift  vector 
is  normp  L  to  the  average  onset-flow  velocity  vector.  The  inlet  and  exit  vel¬ 
ocities  can  be  determined  by  using  (ll)  for  v  and  by  referring  to  the  vector 
diagram  of  figure  3*  All  velocities  will  be  normalized  with  the  modulus  of 
the  average  velocity.  The  inlet  and  exit  velocities  are  written  as  follows: 
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(1?) 


r 

2  s?u 


+  j  £k£_ 
U  J  4SP 


V 

u 


r 

2SPU 


_V  Cl  C 

U  J  ASP 


The  inlet  and  exit  angles  of  attac<  (see  figure  2)  can  be  written  as 


-1 

=  XOJYL 


P  r  -I 

r  L  i 

AIM  OC  +  2  S  P 

=  tojn. 

AWUX+  4SP 

aw.  a 

COd,  OC 

-  1 

OC  =  tam, 
£ 


r 

aim  cl  —  2SP 
COd.  OC 


(13) 


Thus  the  turning  angle  is 


A  =  or-  or  =  Wt 
oc  1  e 


cod.  oc 


/  r  v 

V2SP  / 


CLC 


2SP  GCKbOC 


(14) 


The  average  onset  flow  modulus  U  is  used  in  the  definition  of  all  of  the 
hydrodynamic  coefficients  and  variables  in  the  Douglas  Cascade  Program.  In 
3ome  cases  it  may  be  desirable  to  normalize  the  velocities  involved  with  the 
inlet  velocity  modulus,  Up  In  that  case  we  write 
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The  exit  angle  of  attack  can  be  obtained  from  the  exit  velocity  vector.  Also 
the  turning  angle  can  be  calculated,  since  the  inlet  ar.d  exit  angle  of  attack 
are  known.  To  convert  from  the  system  normalized  with  the  average  velocity  to 
the  one  using  the  inlet  velocity  as  t-e  normalizing  factor,  the  following  set 
of  conversions  can  be  used: 
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7.0  EXAMPLES  AM)  COMPARISONS 


In  order  to  give  an  idea  of  trie  wide  class  of  problems  the  Program  can 
handle,  several  cases  have  been  calculated.  Figures  3>  5,  o,  and  7  present 

a  range  of  conf i rurations  and  geometries  that  can  be  handled  with  ease. 

lb  show  the  accuracy  of  the  Program,  comparisons  with  exact  solutions  and 
experimental  data  are  presented.  Figures  8,  9,  10,  and  11  give  these  com¬ 
parisons  . 

Figure  12  compares  the  Douglas  Cascade  Program  with  a  theory  developed  by 
Garrick,  reference  7- 

These  figures  are  now  described  in  detail. 

7.1  TANDEM  CIRCLES 

To  illustrate  the  fact  that  the  Cascade  Program  can  be  used  for  problems 
not  necessarily  associated  with  the  usual  lifting  cascade,  tne  flow  about  an 
infinite  number  of  nonlifting  bodies  in  tandem  was  calculated  (see  figure  k). 

In  this  particular  case  the  bodies  are  circles;  however,  any  body  can  be  used. 

This  tandem  arrangement  can  be  recognised  as  simply  a  cascade  at  90°angle 
of  attack  rotated  90°  so  that  the  axis  of  repetition,  the  y-axis,  becomes  the 
x-axis.  Recall  that  any  problem  where  the  flow  pattern  repeats  indefinitely 
in  one  direction  can  be  handled  by  the  Douglas  Program. 

7.2  MULTIPLE-CASCADES 

Figures  5  and  6  are  included  to  illustrate  the  multiple -cascade  capability. 

Shown  in  figure  5  are  three  parallel  cascades.  The  central  cascade  has 

circulation,  while  the  other  two  do  not.  It  would  be  a  mistake  to  say  that 
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because  they  have  no  circulation  the  two  outside  cascades  are  nonliftiDg.  The 
resultant  forces  in  an  interaction  problem  are  not  determined  by  the  circula¬ 
tion  alone. 

It  is  true  that  the  total  lift  of  the  entire  set  is  proportional  to  the 
total  circulation  of  the  set  and  that  the  net  drag  of  the  set  is  zero.  However , 
this  is  a  gross  effect  for  the  set  and  does  not  hold  for  the  individual  members. 

As  is  noted  in  figure  5  the  total  drag  coefficient  of  the  first  and  last  cascade 
is  the  negative  of  the  drag  coefficient  on  tne  second  cascade.  The  total  lift 
coefficient  of  the  cascade  set  is  22.7. 

The  trailing  edge  of  the  central  circle  is  at  -30°  from  the  horizontal  and 
the  inlet  angle  of  attack  for  the  cascade  set  is  1*3°. 

For  multiple  cascades  several  limitations  must  be  kept  in  mind.  First,  the 
cascades  must  be  parallel.  Second,  there  can  be  only  one  spacing  associated  with 
all  of  the  cascades.  In  general,  all  the  cascades  of  a  set  must  have  the  same 
spacing;  however,  certain  exceptions  to  this  can  be  made.  One  of  these  exceptions 
is  illustrated  in  figure  6.  In  this  case  the  spacing  of  the  second  cascade  is 
exactly  half  the  spacing  of  the  first.  This  effect  was  obtained  by  putting  two 
.ascade  bodies  in  the  second  cascade.  The  two  bodies  of  the  second  cascade  are 
placed  one  above  the  other  and  spaced  at  exactly  one-half  the  spacing  of  the 
first  cascade.  This  process  can  easily  be  generalized  to  many  cascades,  and  the 
result  is  that  the  spacing  ratio  of  two  or  more  cascades  will  be  a  rational  number. 
It  is  noted  that  there  is  still  only  one  spacing  associated  with  both  cascades 
of  figure  6.  This  spacing  is  shown  in  the  figure. 

Since  neither  cascade  has  circulation,  the  set  has  no  net  lift.  However, 

bodies  2  and  3>  of  the  second  cascade,  have  lift  equal  but  opposite  in  sign  thus 

the  total  lift  is  zero.  The  drag  of  the  first  cascade  is  equal  to  the  thrust  of 

the  second.  The  drag  coefficient  is  0.5^6. 
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7.3  ANALYTIC  TEST  CASES 


The  Douglas  Neumann  Program  with  its  cascade  modification  can  calculate 
the  flow  about  any  cascade  profile.  To  test  this  claim,  the  flow  was  calcu¬ 
lated  about  the  blade  section,  profile  "A",  shown  in  figure  7*  This  extreme 
shape  was  generated  from  a  circle  by  using  a  series  of  conformal  transformations. 
Appendix  B  gives  the  details  of  these  transformations.  The  pressure  coefficients 
obtained  from  the  transformation  method  are  exact.  Figures  8a  and  8b  give  the 
exact  and  calculated  pressure  distributions  over  the  blade  in  a  cascade  at  two 
different  angles  of  attack.  Generally  the  agreement  is  good  for  such  an  ex¬ 
treme  blade  shape;  however,  if  greater  accuracy  were  desired,  more  coordinate 
points  describing  the  body  would  be  needed.  To  prevent  crowding,  not  all  of 
the  points  calculated  by  the  program  are  shown  in  these  figures.  For  an  ex¬ 
ample  where  all  of  the  points  are  plotted  see  figures  12a  and  12b. 

The  blade  section  shown,  profile  "C"  in  figure  9,  was  obtained  through 
a  series  of  conformal  transformations  in  the  same  manner  as  the  blade  of  figure 
7  (see  Appendix  B).  Also,  figures  9a  and  9b  show  the  exact  and  calculated 
pressure  distributions  over  the  blade  in  cascade  at  two  angles  of  attack.  Agree¬ 
ment  between  analytic  and  calculated  pressure  distributions  is  better  than  that 
of  figures  8a  and  8b  as  is  to  be  expected,  since  the  shape  is  less  extreme. 

7.4  EXPERIMENTAL  COMPARISONS 

Shown  in  figures  10a,  10b  and  10c  are  calculated  and  experimentally  ob¬ 
tained  pressure  distributions  for  an  NACA  65-010  cascade  blade  at  three  values 
of  lift  coefficient.  The  experimental  data  were  taken  from  reference  6.  Figure 
11  shows  the  experimental  and  calculated  lift  coefficients  as  functions  of  the 
"effective”  angle-of -attack  for  the  cascade. 

The  "effective"  angle  of  attack  is  simply  the  stagger  angle  plus  the  inlet 


angle  of  attack. 
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The  calculated  and  experimental  pressure  distributiom  agree  quite  well 
except  for  a  small  region  near  the  trailing  edge.  The  aescrepancies  near  the 
trailing  edge  are  probably  due  to  boundary-layer-thicxness  effects. 

7.5  COMPARISON  OF  THE  DOUGLAS  METHOD  WITH  A  METHOD  DUE  TO  I.  E.  GARRICK 

I.  E.  Garrick  (reference  7)  applied  a  straight-line  cascade  transforma¬ 
tion  in  series  with  a  Theodorsen-type  transformation  to  map  an  airfoil  in 
cascade  onto  a  single  circle.  Once  the  transformation  was  obtained,  the  pres¬ 
sure  distribution  over  the  cascade  airfoil  could  be  found.  Accuracy  can  easily 
be  lost  in  the  process  of  carrying  out  the  operations  involved,  because  of  the 
peculiar  nature  of  the  straight-line  cascade  transformation.  This  may  explain 
some  of  the  discrepancies  between  the  two  methods. 

Garrick  calculated  the  flow  over  an  NACA  4412  airfoil  in  cascade.  Shown 
in  figures  12a  and  12b  are  pressure  distributions  calculated  by  Garrick  and 
Douglas  for  the  4412  cascade  at  two  lift  coefficients.  In  figures  12a  and  12b 
all  points  at  which  the  Program  executed  a  calculation  are  shown.  This  will 
serve  as  a  indication  of  the  number  of  coordinates  used  by  tne  Program  in  the 
solution  of  the  problem. 
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APPENDIX  A 


Basic  Solutions 

Eacr  "basic  solution"  is  a  solution  of  the  potential-fluid-flov  problem 
for  a  cascade  of  bodies  with  a  given  onset  flow.  Once  determined,  these  "cas^c 
solutions"  are  combined  in  such  a  way  that  the  desired  flow  at  infinity  is  ob¬ 
tained  and  the  Kutta  condition  on  each  cascade  is  satisfied.  The  "basic  solu¬ 
tions"  needed  for  this  combination  procedure  are: 

1)  Flow  about  tne  several  cascades  at  zero  angle  of  attack, 

2)  Flow  about  the  several  cascades  at  90°  angle  of  attack, 

3)  Flow  about  tne  cascades  due  to  the  presence  of  circulation  in 
the  first  cascade, 

b)  Flow  about  the  cascades  due  to  tne  presence  of  circulation  in 
tre  second  cascade,  and  so  forth. 

To  obtain  circulation  about  a  cascade  profile,  a  unit  vortex  is  placed 
within  tne  profile.  This  vortex  serves  as  the  onset  flow  for  the  circulatory 
"basic  solution". 

In  a  cascade  the  flow  pattern  repeats  indefinitely  along  one  axis,  in  tnis 
case  the  y-axis.  Therefore  it  is  only  necessary  to  deal  with  one  of  the  cascade 
blades.  If  the  Kutta  condition  holds  for  one  blade  of  a  cascade,  it  holds  for 
all  of  them.  Thus,  in  dealing  with  the  cascade,  only  one  blade  will  be  con¬ 
sidered. 

Each  basic  solution  violates  the  Kutta  condition.  A  measure  of  this  vio¬ 
lation  Is  the  difference  AV  of  the  tangential  velocities  above  and  below  the 
trailing  edge.  The  AV* s  of  the  basic  solutions  are  denoted:  AVO,  AV90,  and 
AVmn;  where  m  =  1,  2...NC  and  n  »  1,  2“*NC.  Here  AV^  is  the 
effect  on  the  Kutta  condition  of  the  mth  cascade  due  to  circulation  in  the  nth 
cascade.  NC  is  the  number  of  cascades.  Added  together,  the  basic  solutions 
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must  satisfy  the  Kutta  coniitior  on  eacr.  cascade  and  also  give  the  desired 


flow  at  infinity. 


COMBINATION  OF  BASIC  SOLUTIONS 


For  a  set  of  isolated  airfoixs  the  following  set  of  linear  equations  in 
the  unknown  circulation  strengths  satisfies  the  Kutta  conditions  on  each  air¬ 
foil.  The  uniform  onset  flow  at  infinity  is  of  speed  U  and  angle  of  attack  a. 
The  set  of  linear  equations  with  the  unknown  rm  is 

AVO-jUcosa  +  AV90^Usina  +  AV^  +  AV^2  P  2  +  •**  =  0 

A  VC^Ucosa  +  AV^OjUsina  '  AV^  Hl  +  AV22  ^2  +  *  *  *  11  0 

*  • 

•  • 

#  • 

AVONCUcosa  +  AV90NCUsina  +  AVNC1  ^1  +  AVNC2  *~2  +  ***  =  0 


or,  in  matrix  form. 


■ 

'  > 

Ayn  AV12  •  •  • 

rl 

AVO^cosa  +  AV90}Sina 

AV?1 

r2 

AVOpCosa  +  AV902sina 

• 

• 

< 

• 

* 

►  =  -U  < 

> 

• 

* 

A%C1 

^NC 

AVONCcosa  +  AV90NCsina 

- 

- 

(Al) 


where 

a  =  angle  of  attack  of  the  set  of  airfoils 
P  =  circulation  strength  per  airfoil 

For  a  set  of  cascades  there  is  a  complicating  factor,  and  (Al)  cannot  be  used 
directly.  The  angle  a,  the  average  angle  of  attack,  is  not  necessarily  known. 
V.'hat  is  known  is  any  one  of  the  following: 

1)  Cascade  or  set  of  cascades  at  a  prescribed  average  angle  of  attack  (a). 

2)  Cascade  or  set  of  cascades  at  a  prescribed  inlet  angle  of  attack  (atj). 

^ )  Cascade  or  set  of  cascades  at  a  prescribed  li  ft  coefficient  (Cp,). 

4)  Cascade  or  set  of  cascades  wit-  a  prescribed  turning  angle  (Aa). 

22 


These  cases  are  mutually  exclusive.  For  example:  if  Oj-  is  prescribed  the 
other  three  cannot  be  prescribed. 


For  case  (l)  equation  (Al)  can  be  used  directly.  For  cases  (2),  (3),  and 
(4)  there  is  an  additional  unknown,  namely,  a.  Thus  an  additional  equation 
must  be  found  for  these  cases.  The  additional  equations  needed  for  cases  (2), 
(3),  and  (4),  respectively,  are:  MC 


OKI  OC  'tdJAj  0CT  -  A aJYU  OC  -  -J—  > 

1  2SP  /—j 


p  =  0 
m 


wc 

Yr.-Q.o 


m-  1 


2. 

Uc  /  'm 
m=  1 


wc 


C£ra 


~tcun.  (A  Oc)  - 


SP 


rP- 


m  =. » 
THT 


(A2-2) 


(A2-3) 


( A2-4) 


j_Y  r 

ZSP  Z-  m 
m= 1 

With  some  rearranging  of  the  linear  equations,  (A 2-2)  and  (A2-3)  can  be  in¬ 
corporated  into  the  set  of  linear  equations  (Al).  However,  (A2-4)  cannot  be 
incorporated  and  must  be  solved  iteratively  together  with  (Al). 


When  oij  is  the  desired  input,  equation  (A2-2)  can  be  incorporated  into 
(Al).  The  resulting  matrix  equation  is 
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When  the  lift  coefficient  is  input,  equation  (2-3)  can  be  used  wit;, 
(Al)  in  the  following  manner.  First  we  may  write 
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UAV901  AV1]l 


AV1?  *  '  '  AV1NC  i  f  tana 


UAV9°NC  AVNC1  A^ 


'NCNC 


P  ^/cosa 
P  j/cosa 


r  Nc/i 


cosa 


-UAVO-, 


-UAVO, 


-IAVOjic 


'cosa 


-(AH) 


In  the  above  set  of  linear  equations  notice  tnat  cosa,  with  a  an  ur.Known, 
still  appears  on  the  r.pht-hand  side  of  tie  equation.  To  solve  this  set  of 
equations,  define  A  as 

f-V90i  AVU  AV  •  .  -"I 


A  = 


UAV^c  avnc1  *  * 


Let  A  ,  the  inverse  of  A,  be  defined  as 
f  all  ai2  *  ’  ’  alNC  "1 


Then  (M)  becomes 


tana 


aNCNC 


-UAVO, 


'cosa 


-UAVOl 


P  NC/ 


cos a 


and  therefore 


-UAVOnc 

Cl/cosa 


tana  “  AVOi  an  +  AVO?  a12  +  *  *  *(CL/cosa^  a] 
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Equation  ( A5 5  can  then  he  used  to  solve  for  cosa.  Equations  (Ah)  and  (A5) 
represent  the  solution  vhen  is  input. 

As  has  been  stated  before,  (Al)  and  (A2-^)  must  be  solved  iteratively 
vhen  the  turning  angle  £a  is  desired. 

When  dealing  vith  the  usual  case  of  a  single  cascade,  the  equations 
become  very  simple.  Equation  (Al)  reduces  to  one  equation  ani  the  matrix  of 
(A3)  is  of  order  two  and  its  solution  is  trivial. 
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APPENDIX  B 

Analytic  Cascade  Test  Cases 

In  order  to  check  the  Douglas  Cascade  Program,  several  exact  cascade 
solutions  were  generated  by  conformal  transformation  methods.  Specifically, 
a  circle  was  mapped  by  a  series  of  transformation  functions  into  a  profile 
in  cascade.  The  final  shape  could  not  be  determined  exactly  ariead  of  time 
but  certain  characteristics  could  be  controlled. 

The  first  in  the  series  of  mapping  functions  was  the  Karman-Trefftz 
transformation.  Tnis  function  maps  a  circle  in  the  complex  S-plane  onto  an 
airfoil  shape  in  the  complex  Q-plane.  The  transformation  is 


The  final  airfoil  shape  is  determined  by  the  location  of  the  circle  with 
reference  to  the  S  coordinate  system.  The  real  constant  d  is  the  distance 
from  the  origin  to  the  intersection  of  the  circle  with  the  real  S 
axis.  This  intersection  maps  to  the  trailing  edge  of  the  airfoil.  This  mapping 
is  very  similar  to  the  Joukowski  transformation,  in  that  a  displacement  of  the 
circle  toward  the  negative  real  S-axis  produces  thickness  and  a  displacement 
toward  the  positive  imaginary  S-axis  produces  positive  camber.  Thus  tnese 
three  parameters  determine  the  Karman-Trefftz  airfoil:  (l)  the  trailing-edge- 
angle  constant,  r;  (2)  the  radius  of  the  circle,  a;  and  (3)  the  location  of 
the  center  of  the  circle  in  the  S  coordinate  system.  For  profile  "A"  r  *  1.8, 


26 


a  =  1.0  and  the  center  location  is  (-0.1,  1  /  ^j~2  i).  For  profile  "C" 
r  =  1.85,  a  =  1.0,  and  the  center  locaton  is  {-0.02,  0.5i). 


The  second  and  last  transformation  takes  the  Karman-Trefftz  airfoil  in 
the  Q-plane  into  some  profile  in  a  cascade  in  the  final  z-plane  by  using  tne 
following  transformation: 


z 


(B3) 


The  spacing  of  the  resulting  cascade  is  always  2«,  and  A  and  B  are  complex 
quantities.  The  singular  points  at  A  and  B  must  be  outside  of  the  airfoil  and 
must  not  touch  its  surface.  Essentially,  the  points  A  and  B  in  the  ^-plane 
go  to  plus  and  minus  infinity,  respectively,  in  the  z-plane.  The  closer  the 
points  are  to  the  airfoil  surface,  the  longer  the  cascade  profile.  For  con¬ 
venience  of  computation,  the  points  A  and  B  are  not  selected  in  the  airfoil 
plane  but  in  the  circle  plane.  In  the  circle  plane  they  are  called  A'  and  B’ . 
For  profile  "A" 

A1  =  0.c57  +  0i,  B'  =  -1.1  +  O.-i.  For  profile  "C" 

A'  =  1.0 j60-  +  Oi ,  B'  =  -1.08to0-;  +  Oi .  To  determine  the  coordinates 

of  A  and  B  the  transformation  of  (B6)  is  used  with  A*  and  B'  in  place  of  S 
and  A  and  B  in  place  of  Q. 

To  obtain  t..e  flow  field  in  tne  z-plane,  tne  complex  potential  must  be  differ¬ 
entiated  . 

If  F  is  this  complex  potential  and  w  the  complex  velocity, 

j,)  .  ,BU) 

'  dz  dS  dQ  dz 

dS  .  , 

The  term  cTQ  can  be  obtained  from  (Bl) 


o  O 

dS  =  Sc  -  d^ 
dQ  q2  _  r2d2 


(B5) 
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The  derivative 


oFiS} 


is  the  complex  velocity  in  the  circle  plane;  i.e., 


it  represents  the  flov  field  about  the  circle.  The  question  to  consider  is: 


in  what  flow  field  is  the  circle  immersed?  The  answer  lies  in  the  transformation 


(B3),  '  ,e.,  tie  mapping  function  that  takes  tne  Karman-Trefftz  airfoil  to  a  pro¬ 
file  in  cascade.  It  can  be  shown  tnat  a  source  at  B  and  a  sink  at  A  in  the 
airfoil  plane  r've  a  uniform  flow  from  minus  to  plus  infinity  in  the  cascade 
plane.  It  can  also  be  shown  tnat  a  vortex  at  each  of  tnese  points  gives  a 
vertical  component  to  tne  uniform  flow  ir  the  cascade  plane.  Therefore  the 
following  relations  hold: 

U-jCosccT  =  source  strengt-  ,  v,  at  B 

Lisinas  =  vortex  strer.ctn,  P ,  at  B  ( B  ’ ) 

U,  cosa.  =  sink  strength,  M,  at  A 

•UgSina^  -  vortex  strength,  T  ,  at  A 


To  preserve  continuity  of  miss  in  the  cascade  plane,  tie  source  strength 
at  B  mus  t  be  equal  to  v  e  s:r.K  strength  at  A.  Tne  flow  in  which  tne  cyl.nuer 
is  immersed  is  trier  -enerated  by  a  source  and  vortex  at  B  and  a  sink  and 
vortex  at  A. 


It  can  be  s' own  that  the  circulation  about  tne  cascade  profile  is  tie  r.e  -- 
ative  difference  of  the  strenptns  of  the  vortexices  located  a.  A  and  B.  T  us 
all  irculations  add  uh  to  zero.  The  onset  flow  :,as  now  beer;  determined  and 
therefore  the  flow  about  a  circle  in  this  cnset  flow  can  be  determined.  Tne 
complex  potential  function  for  this  flow  is 


J 


F  =  K,  ln(S-A)  +  K-i  Inf  S 


-  =-) 
A 


Kg  ln(S-B)  +  ^  Ln(S  -  |-)  (b8) 


where  a  is  the  circle  radius  and  K  is  the  complex  source  strength  M  +["1. 
M  is  the  source  mass  flow  and  T  is  the  circulation.  The  complex  velocity  is 
just  the  derivative  of  F  and  is 


dF 

dS 


(B9) 


Ki  =  U-j^cosOj.  +  i  sinaj),  K. 


U 

E  , 

-Uj(cosctj  +  i  y-  sinaE) 


If  t:e  expressions  for  the  deri”ities  of  (B  ),  (Br)  and  (B9)  are  substi¬ 
tuted  into  equation  (A9)  and  if  S  ta-.es  the  values  of  the  coordinates  of  the 
circle  tr.e  velocity  over  tre  surface  of  trie  cascade  body  in  the  z -plane  can  be 
determined.  In  tr.ese  formulas  tr.e  velocities  are  normalized  with  the  inlet 
velocity  modulus  U-. 


29 


APPENDIX  C 


Program  Input  and  Output 

A  summary  of  the  Program  input  is  presented  before  the  detailed  expiate 
tion  Is  given. 

As  is  stated  in  the  text,  one  or  more  cascade  bodies  of  arbitrary  shape 
can  be  handled  by  the  Program.  However,  there  are  three  practical  restrict¬ 
ions  on  the  input,  as  follows: 

vl)  The  cascade  body  or  bodies  must  be  of  finite  thickness. 

(2)  The  maximum  number  of  bodies  with  circulation  is  8. 

(3)  The  maximum  number  of  points  describing  all  of  the  bodies  and 
off-body  points  is  500. 

Each  body  considered  is  a  lifting  body  with  a  stagnation  point  at  the 
trailing  edge,  unless  otherwise  specified  in  the  input.  It  Is  assumed  in  the 
program  that  the  first  coordinate  point  input  is  the  trailing  edge.  The  sur¬ 
face  coordinates  are  input,  starting  from  the  trailing  edge  and  progressing 
around  the  body  in  the  clockwise  direction.  The  last  point  of  a  body  must  be 
the  first  point  repeated.  However,  if  the  body  is  ron-closed  and  non-lifting 
the  last  point  i3  not  the  first  point  repeated. 

In  addition  to  calculating  the  flow  on  the  cascade-body  surface,  the 
Program  can  calculate  the  flow  at  points  in  the  flow  field.  The  coordinates 
of  these  off -body  points  are  input  in  the  same  manner  as  the  coordinates  of 
the  cascade  body. 

All  coordinates,  whether  bod„ -surface  or  off-body  coordinates,  may  be 
scaled,  rotated,  and  translated.  The  Program  executes  these  operations  in 
the  order  named. 
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For  diagnostic  purposes  the  two  matrices  A 


can  be  printed  out. 


Input  Data 


and 


Each  case  must  consist  of  a  header  card,  case  control  data,  body  control 
data,  and  coordinate  data.  The  header  card  contains  the  description  of  the 
case,  control  flags,  and  case  number.  The  case  control  data  specify  certain 
constants  used  in  the  computation.  The  body  control  data  specify  the  amount 
of  coordinate  data  being  input  and  constants  used  to  modify  the  coordinate 
data.  The  coordinate  data  describe  either  the  tvo-dimensional  cross  section 
of  the  body  or  off -body  points.  All  data  cards  must  contain  sequence  numbers 
in  card  columns  77  through  80,  so  that  the  data  may  be  sorted.  If  any  data 
cards  are  found  to  be  out  of  sequence,  the  program  will  discontinue  execution. 

The  data  must  be  arranged  in  the  following  order  (see  figure  2): 

1.  Header  Card  (l  card) 

2.  Case  Control  Card  (I  card) 

3.  Body  Control  Cards  (2  cards) 

k.  Coordinate  Data  Cards  (Variable  number  of  cards) 

Items  3  and  4  are  repeated  for  each  body  (if  more  than  one  Jody  is  being 
considered)  and  for  off -body  points,  if  any.  The  Y  coordinates  for  each 
body,  or  off -body  points,  must  always  start  on  a  new  card,  and  must  always 
follow  the  X  coordinates.  Additional  cases  may  be  run  by  placing  additional 
sets  (items  1  through  4)  one  after  another. 


If  additional  cases  are  to  be  run  using  some  or  all  of  the  previous 
untransformed  coordinate  data,  the  "Subcase”  capability  is  used.  All  that 
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need  be  Input  are  the  body  control  cards,  vith  the  "Subcase”  flag  marked, 
that  will  transform  the  coordinates  of  the  previous  case  in  a  manner  desired 
for  the  present  case. 

If  additional  bodies  are  to  be  added  their  coordinates  are  input  in  the 
normal  manner,  without  the  "Subcase"  flag,  following  the  body  control  cards 
mentioned  above.  If,  in  additional  cases,  bodies  are  to  be  deleted  or  re¬ 
placed  these  bodies  must  appear  last  in  the  sequence  of  bodies  in  the  original 
case.  To  delete  body  coordinates  simply  omit  the  body  control  cards  for  that 
body. 

When  replacing  a  body  for  a  subcase  simply  introduce  the  new  body  con¬ 
trol  cards  and  coordinates  in  the  place  of  the  replaced  body. 

As  an  illustration,  if  it  is  desired  to  run  two  Cascades,  call  them  A 
and  B,  and  then  to  delete  B  and  run  A  alone,  or  with  a  new  body  C,  the  following 
procedure  is  followed. 

Input  A  as  the  first  body  and  B  as  the  second  as  shown  in  figure  2.  Then 
write  new  header  and  case  control  cards  for  the  second  case, placing  them  in 
back  of  the  y  cards  for  body  B.  The  body  control  cards  for  A  are  written 
with  the  "Subcase"  flag  marked.  To  omit  body  B  simply  omit  the  body  control 
cards  for  body  B.  If  a  new  body  C  is  to  be  input  in  place  of  B, its  body  con¬ 
trol  cards  followed  by  its  coordinate  cards  would  be  input  following  the  body 
control  cards  for  A. 

A  second  type  of  subcase  capability  exists.  If  only  the  case  control 
data  is  to  be  changed  for  a  second  case  simply  mark  the  flag  in  card  column  8 
and  write  out  the  new  case  control  data.  For  example,  if  the  calculation  is 
desired  at  a  second  value  of  CL  simply  write  an  additional  header  card  and 
case  control  card  with  the  new  value  of  CT  .  This  may  be  repeated  indefinitely. 
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~y5\  /  “  Y  Coordinates 

/  -  X  Coordinates 
-  Body  Control  Card  2 


-  Body  Control  Card  1 


-  Y  Coordinates 


[/  -  X  Coordinates 
-  Body  Control  Card  2 
-  Body  Control  Card  1 
-  Case  Control  Card 
~  Header  Card 


Figure  2.  -  Arrangement  of  data  cards 


Every  complete  job  must  be  followed  by  an  end-of -file-card  (7-8  punch 
in  card  column  l)  and  a  ^IBSYS  card  (^IBSYS  in  card  columns  1-6).  Another 
case  may  be  run  by  placing  a  header  card,  case  control  card,  etc.,  after  the 
last  Y  coordinates  of  the  first  case.  Although  not  part  of  the  data,  an 
end-of-file-card  must  always  follow  the  program  deck,  so  that  the  data  for 
any  job  are  always  preceded  by  and  followed  by  end-of-file  cards. 

Header  Card: 

Card  column  1  must  always  be  filled  with  any  nonzero  integer  that  indi¬ 
cates  the  number  of  bodies  being  input.  This  integer  must  not  be  larger  than  8. 
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Card  columns  2-12,  when  punched  with  any  nonzero  integer,  activate  flags  that 
indicate  the  following: 


Cerd  column: 
2 

3 

k 

5 


6 

7 

8 


9-12 

13-60 


63-68 


77-80 


Flow  is  to  he  determined  at  points  off  the  "body, 
a  is  input  for  use  in  the  combination  equations  for  airfoils. 
Act  is  input  for  use  in  the  combination  equations  for  airfoils. 
Inlet  a  is  input  for  use  in  the  combination  equations  for 
airfoils . 

C_  is  input  for  use  in  the  combination  equations  for  airfoils. 
The  matrix  of  influence  coefficients  is  to  be  printed  out. 

Go  directly  to  combination  solution  using  basic  velocity  solu¬ 
tions  of  the  previous  case. 

Hot  used. 

This  description  of  the  case  will  be  printed  on  each  section 
of  output. 

This  case  number  will  be  printed  at  the  beginning  of  the 
output . 

A  sequence  number  must  appear  in  these  columns. 


Case  Control  Data: 

All  of  ths  input  items  defined  in  this  section,  with  the  exception  of 

CHORD,  are  assumed  to  be  zero  if  no  value  is  input. 

CHORD  The  chord  length  to  be  used  in  computations  for  this  case.  It 

will  be  assumed  to  be  1.0  if  no  value  is  input.  Any  value 
input  must  appear  with  a  decimal  point. 

SPACING  The  spacing  between  the  bodies  of  the  infinite  cascade.  A 

decimal  point  must  be  specified. 


3^ 


a 

INLET  a 

Ha 


Body  Control  Data 
Body  Control  Card 
NN 


MX 

MY 

THETA  -  Stagger 
Angle 


The  lift  coefficient  to  be  used  in  the  combination  equations 
for  airfoils.  A  decimal  point  must  be  specified. 

The  angle  of  attack  (in  degrees)  to  be  used  in  the  combina¬ 
tion  equations  for  airfoils.  A  decimal  point  must  be  speci¬ 
fied. 

The  inlet  angle  (in  degrees)  to  be  used  in  the  combination 
equations  for  airfoils.  A  decimal  point  must  be  specified. 
The  cascade  turning  angle  (in  degrees)  to  be  used  in  the 
combination  equations  for  cascades.  A  decimal  point  must 
be  specified. 


1 

The  number  of  points  on  the  body  being  input.  The  sum  of 
all  the  NN  for  all  bodies  in  a  case  must  not  be  greater 
than  ‘,>00.  This  number  must  not  specify  a  decimal  point, 
and  it  must  be  punched  at  the  far  right  of  its  field  (right- 
justified). 

The  factor  used  to  multiply  all  x  coordinates.  It  is 
assumed  to  be  1.0  if  no  value  is  input.  A  decimal  point 
must  be  specified. 

The  factor  used  to  multiply  all  y  coordinates.  Otherwise, 
same  as  MX . 

The  angle  (in  degrees)  through  which  all  points  are  to  be 
rotated  about  the  origin  in  tne  clockwise  direction  (stagger 
angle).  A  decimal  point  must  be  specified. 


35 


ADDX  The  constant  to  be  added  to  all  x  coordinates.  A  decimal 

point  must  be  specified. 

ADDY  Tne  constant  to  be  added  to  all  y  coordinates.  A  decimal 

point  must  be  specified. 

Body  Control  Card  2 

BDH  The  body  sequence  number.  This  number  must  be  a  nonzero 

integer  if  body  coordinates  follow;  it  is  zero  only  if  off- 
body  coordinates  follow. 

HU  This  is  a  flag  that  must  be  any  nonzero  integer  only  if  the 

body  whose  coordinates  follow  is  to  be  considered  a  nonlift¬ 
ing  or  noncirculatory  body. 

SUBCASE  This  is  a  flag  that  directs  the  program  to  use  the  unmodified 

coordinates  of  the  body  of  the  previous  case.  It  must  be  any 
nonzero  integer. 

XMC  The  coordinates  of  the  moment  center  to  be  used  vhen  the 

ym: 

combination  equations  for  airfoils  are  used.  A  decimal 
point  must  be  specified. 

Coordinate  Data: 

In  inputting  the  body  coordinates,  it  is  essential  that  the  coordinate  data 
start  at  the  body  trailing  edge,  progress  around  the  body  in  the  clockwise 
direction  and  that  the  last  point  input  be  the  first  point  repeated  for  a 
closed  body  only.  The  example  problem  illustrates  this  procedure  in  figure 
13b  and  13c. 

X  The  x  coordinates  of  the  points  defining  the  body,  (The 

x  coordinates  of  the  off -body  points  if  BDN  is  zero).  A 
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decimal  point  must  be  specified. 

The  y  coordinates  of  the  points  defining  the  body.  (The 
y  coordinates  of  the  off -body  points  if  BDH  is  zero).  A 
decimal  point  must  be  specified. 

The  x  coordinates  must  precede  the  y  coordinates  in  the 
deck  arrangement.  There  must  be  US  x  coordinates  and  NN 
y  coordinates . 

Output ; 

All  sections  of  output,  vith  the  exception  of  the  matrix  printout,  is 
preceded  by  the  following  header: 

DOUGLAS  AIRCRAFT  COMPANY 
LONG  BEACH  DIVISION 

The  information  contained  on  the  header  card  and  the  case  control  card  is 
printed  on  the  first  page  of  output.  On  this  page,  FLAG  2  corresponds  to 
card  column  2  on  the  header  card,  FLAG  3  bo  card  column  3,  and  so  on. 

The  next  section  of  output  consists  of  basic  data,  i.e.,  control  and 
coordinate  data,  for  each  body.  The  body  control  data  are  printed  out  first 
the  column  headers  follow: 

X  Y  DELTA  S  SUMDS  D  ALPHA 

The  X  and  Y  columns  list  the  modified  x  and  y  coordinates  and 
the  midpoints  of  the  element  formed  by  two  consecutive  body  points.  A  modi¬ 
fied  coordinate  is  one  that  has  been  scaled,  ”otated  and  translated  accord¬ 
ing  to  the  input.  The  column  headed  by  DELTA  S  specified  the  length  of  the 
element  formed  by  two  consecutive  body  points,  SUMDS  shows  a  running  sum 

of  the  DELTA  S  column,  and  D  ALPHA  shows  the  angle  between  two  consecu¬ 
tive  elements. 
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If  off-body  points  are  input,  the  following  column  headers  are  printed 
out  after  the  basic  data  for  all  of  the  bodies: 

X-OFF  Y-OFF 

These  columns  merely  list  the  modified  off -body  points. 


If  FLAG  7  i8  punched  (card  column  7  on  the  header  card),  the  A  and 

JK 

B  matrices  are  printed  cut  after  the  basic  data.  Both  matrices  are  printed 

jK 

out  in  row  order  across  the  page. 


If  off -body  points  are  being  run,  A.,  and  B  off-body  matrices  are 

JK  JK 

formed  and  are  printed  out  after  the  on-body  matrices. 


The  next  section  of  output  consists  of  the  original  x  and  y  coordin¬ 
ates  and  the  midpoints  of  the  elements  formed  by  these  coordinates,  the 
velocities  at  the  midpoints  of  the  elements  (V),  and  the  corresponding  pres¬ 
sure  coefficients  (CP)  and  source  densities  (SIGMA).  This  output  is  repeated 
for  each  onset  flow  i.e.,  basic  solution. 


The  combination  solution  follows  the  basic  so.lutions.  The  combination 
solution  is  a  suitable  combination  of  the  basic  solutions  in  such  a  way  that 
the  Kutta  condition  is  met  and  the  other  input  requirements  satisfied.  The 
format  is  exactly  the  same  as  that  of  the  preceding  solutions  for  one  onset 
flow,  except  that  the  SIGMA  column  is  replaced  by  a  DELTA  S  column,  which 
specifies  the  lengths  of  the  original,  unmodified  elements.  Also  shown  are 
computed  and  input  constants  that  apply  to  the  combination  solution:  spacing, 
alpha,  X^,  Y  ,  inlet  alpha,  exit  alpha,  inlet  velocity,  exit  velocity, 
and  delta  alpha. 


The  last  section  of  output  shows  the  solution  at  the  points  off  the 
body.  The  column  header  for  this  section  is  simply: 
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XT  VXL  VYL 

where  the  X  and  Y  columns  list  the  off -body  points  and  the  VXL  and 
VYL  columns  list  the  X  and  Y  components  of  velocity  at  the  specified 
off -body  point. 

Example  Problem 

To  illustrate  ohe  input  procedure  and  to  shov  and  example  of  computed 
output,  an  example  problem  is  presented.  The  problem  consists  of  a  cascade 
of  circles  of  unit  radius,  spaced  three  radii  apart.  The  cascade  is  at  an 
average  angle  of  attack  of  10°.  The  only  coordinate  modification  is  a  rota¬ 
tion  of  180°  to  place  the  first  coordinate  point  at  the  desired  trailing- 
edge  position  which,  in  this  problem,  is  on  the  x-axis.  (see  figure  13a) 

The  circle  is  composed  of  thirty  coordinate  points  spaced  equally  around  the 
perimeter.  The  input  coordinates  (see  figure  13b  and  13c)  progress  in  the 
clockwise  direction. 

The  example  problem  also  shows  the  input  and  output  for  four  off -body 
points.  The  coordinates  of  these  points  are  input  in  the  same  manner  as  the 
circle  coordinates  (see  figures  13d  and  13e)  except  that  the  flag  BDN  is 
marked  0. 

The  output  for  the  example  problem  is  shown  in  figures  l4a  through  f. 
Figure  l4a  shows  the  basic  data,  i.e.,  the  input  data  and  the  transformed 
coordinates.  The  basic  data  are  given  on  the  first  three  pages  of  output. 
Figures  14b,  c,  and  d  show  the  output  sheets  that  give  the  three  basic  flow 
solutions:  the  solution  at  0°  angle  of  attack,  the  solution  at  90°  angle  of 
attack,  and  the  solution  due  to  a  circulatory  flow.  Figure  l4e  shows  the 
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output  sheet  that  gives  the  combination  solution  for  an  average  angle  of 
attack  of  10°.  Also  shown  on  the  output  of  figure  Ike  are  the  following: 

(1)  Inlet  and  exit  velocity  and  angle  of  attack 

(2)  Cascade  lift  and  moment  coefficient  and  the  x  and  y  force 
coefficients 

(3)  Spacing 

(4)  Cascade  turning  angle 

Figure  l4f  shows  the  output  sheet  for  the  off-body  point  velocities. 

Figure  15  gives  a  complete  FORTRAN  IY  listing  of  the  Douglas  Cascade 
Program.  Ten  tape  units  are  needed  on  the  computer  for  this  program. 
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CASCADE  2 


Figure  3.  -  Vector  diagram  showing  inlet,  average,  and  exit  velocities  and  angles  of  attack.  The 
cascade  parameters,  spacing  and  stagger  are  also  shown. 
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Figure  5.  Pressure  distributions  on  three  circles  in  cascade.  The  central  cascade  has  circulation 
SP  =  b  chord,  a,  =  43°,  C,  =  -0.08,  Cn  =  1.36,  C,  =  19.22,  Cn  =  -1.44,  C,  =  3.41 


Figure  6.  Pressure  distributions  on  three  circles  in  cascade,  ’‘he  arrangement  shown  simulates 
two  cascades  with  a  spacing  ratio,  SP,/SP2  =  2.  Both  cascades  are  non-circulatory. 
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Figure  8.  Comparison  of  analytic  and  calculated  pressure  distributions  on  profile  'A"  in  cas¬ 
cade.  (a)  CL  =  0,  6  =  0°,  a,  =  -45°,  SP  0.538 
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Figure  9.-  Continued  (b)  CL  =  0,  6  =  0°,  a,  =  -30°,  SP  =  0.795 
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Figure  10.  -  Comparison  of  calculated  and  experimental  pressure  distributions  on  an  NACA  65-010 
airfoil  in  cascade,  (a)  CL  =  -0.135,  0  -  -33°/  ai  =30°,  SP  =  1.0.  (b)  CL  =  0.2„ 

0  =  -21°,  a,  =  30°,  SP  =1.0 
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Figure  11.  Comparison  of  calculated  and  experimental  lift  coefficient  versus  ‘effective*  angle 
of  attack  for  the  NACA  65-010  airfoil  in  cascade. 
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Figure  12.  Comparison  of  the  pressure  distribution,  as  calculated  by  I.  E.  Garrick  and  by  Douglas 
of  an  NACA  4412  airfoil  in  cascade,  (a)  CL  -  1.0,  8  =  0°,  SP  =  0.968 

(b)  C,  =  1.0,  6  =  45°,  SP  =  1.096 
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Figure  13.-  Example  problem  control  and  data  input  sheets,  (a)  Header  card  and  case  control  daia. 

(b)  Control  data  and  x  coordinates  for  the  cascade  body,  (c)  y  coordinates  for  cascade  body 
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Figure  13.  -  Continued  (d)  Control  data  and  x  coordinates  for  off  body  points. 

(e)  y  coordinates  for  off-body  points. 
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Figure  14.-  Program  output  sheets  for  example  problem,  (a)  Input  or  basic  data. 
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Ic.- 

-j-^cJr  bQLj-tr 


BEST 

AVAILABLE  COPY 


DOUGLAS  AIRCRAFT  COMPANY 
LONG  BEACH  DIVISION 


2"0  CASCADE  TEST  PROBLEM 
CASE  ONE 

STREAMFLOW  SOLUTION 
UNTRANSFORMED  COORDINATES 


I 

X 

-0.99999999 

Y 

-0.00000001 

V 

CP 

SIGMA 

2 

-0.98907378 

-0.97814759 

0.10395583 

0.20791168 

-0.27628037 

Oj.  92366915 

-3.13563100 

3 

-0.94584651 

-0.91354544 

0.30732416 

0.40673663 

-0.82911973 

0.31256047 

-3.05742222 

A 

-0.86128122 

-0.80901700 

0.49726093 

0.58778523 

-1.33354529 

-0.91419758 

-2.89666861 

5 

-0.73907379 

-0.66913059 

0.66546501 

0.74314480 

-1.94094200 

-2.76725584 

-2.63964686 

6 

-0.58456529 

-0.49999999 

0.80458509 

0.86602539 

-2.49449280 

-5.22249430 

-2.25525227 

7 

-0.40450849 

-0.30901700 

0.90854094 

0.95105650 

-3.00859588 

-8.05164909 

-1.69504222 

8 

-0.20677274 
-0. 10452846 

0.97278919 

0.99452189 

-3.39735663 

-10.54203200 

-0.92603669 

9 

0. 

0.10452846 

0.99452189 

0.99452189 

-3.54606095 

-11.57454824 

-0.00000057 

10 

0.20677272 

0.30901699 

0.97278919 

0.95105650 

-3.39735693 

-10.54203403 

0.92603564 

11 

0.40450849 

0.50000000 

0.90854095 

0.86602540 

-3.00859627 

-e. 05165148 

1.69504154 

12 

0.58456530 

0.66913061 

0.80458511 

0.74314483 

-2.49449351 

-5.22249788 

2.25525099 

13 

0. 739073  79 
0.80901697 

0.66546503 

0.58778524 

-1.94094250 

-2.76725775 

2.63964698 

14 

0.86128121 

0.91354544 

0.49726094 

0.40673663 

-1.38354616 

-0.91419996 

2.89666799 

15 

0.94584652 

0.97814760 

0.30732416 
0.20791 169 

-0.82912060 

0.31255904 

3.05742201 

16 

0.98907379 
0. 79999999 

0.10395584 

0. 

-0.27628082 

0.92366891 

3.13563150 

17 

0.98907379 

0.97814760 

-0.10395584 
-0.20791 169 

0.27628004 

0.92366934 

3.13563138 

18 

0.94584651 

0.91354544 

-0.30732416 

-0.40673663 

0.82911979 

0.31256037 

3.05742210 

19 

0.86128121 

0.80901699 

-0.49726094 

-0.58778524 

1.38354519 

-0.91419728 

2.89666864 

20 

0.73907379 

0.66913059 

-0.66546503 

-0.74314483 

1.94094194 

-2.76725560 

2.63964716 

21 

0.58456530 

0.50000000 

-0.8045851 1 
-0.86602540 

2.49449295 

-5.22249502 

2.25525236 

22 

0.40450849 

0. 30901697 

-0.90854095 

-0.95105650 

3.00859606 

-8.05165017 

1.69504198 

23 

0.20677272 

0.  10452846 

-0.97278919 

-0.99452189 

3.39735678 

-10.54203308 

0.92603677 

24 

-0. 

-0.  10452846 

-0.99452189 

-0.99452189 

3. 546061 1 3 

-11.57454944 

0.00000054 

25 

-0.20677273 

-0.30901700 

-0.9727R91 9 
-0.95105650 

3. 39735699 

-10.54203439 

-0.92603590 

26 

-0.40450849 

-0.50000000 

-0.90854094 

-0.86602539 

3.00859636 

-8.05165195 

-1.69504163 

27 

-0.58456-530 

-0.66913061 

-0.80458509 

-0.74314480 

2.49449340 

-5.22249728 

-2.25525141 

28 

-0. 73907380 
•  0.80901700 

-0.66546501 

-0.58778523 

1.94094259 

-2.7672581 1 

-2.63964647 

29 

•0.86128122 

-0.91354544 

-0.49726093 

-0.40673663 

1.38354608 

-0.91419975 

-2.89666826 

30 

-0*94584651 

-0.97814760 

-0.30732416 

-0.20791168 

0.82912064 

0.31255897 

-3.05742183 

31 

-0.98907375 

-0.99999990 

-0. 10395583 
0.00000001 

0.27627993 

0.92366940 

-3.13563240 

Figure  14.  -  Continued  (b)  Solution  at  0°  angle  of  attack 
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DOUGLAS  AIRCRAFT  COMPANY 
LONG  REACH  DIVISION 


2-0  CASCAOE  TEST  PROBLEM 


CASE  ONE 

90-0EGRF.E  FLOW  SOLUTION 
UNTRANSFORMED  COORDINATES 

X  Y  V 


l 

-0.99999999 

-0.000000Q1 

-1.35587278 

-0.98907378 

0.10395583 

-1.53488527 

2 

-0.97814759 

0.20791168 

-1.08526984 

-0.94584651 

0.30732416 

-1.44404635 

3 

-0.91354544 

0.40673663 

-0.61303900 

-0.86128122 

0.49726093 

-1.27005473 

4 

-0.80901700 

0.56778523 

-0.06059517 

-0. T3907379 

0.66546501 

-1.02985202 

5 

-0.66913059 

0.74314480 

0.43530724 

-0.58456529 

0.80458509 

-0.75146042 

6 

-0.49999999 

0.86602539 

0.77778012 

-0.40450849 

0.90854094 

-0.47140203 

7 

-0.30901700 

0.95105650 

0.95146379 

-0.20677274 

0.97278919 

-0.22030935 

a 

-0. 10452846 

0.99452189 

0.99999999 

0. 

0.99452189 

-0.00000036 

9 

0.10452846 

0.99452189 

0.95146409 

0.20677272 

0.97278919 

0.22030865 

10 

0.30901699 

0.95105650 

0.77778073 

0.40450849 

0.90854095 

0.47140139 

ll 

0.50000000 

0.86602540 

0.43530802 

0.58456530 

0.80458511 

0.75145990 

12 

0.66913061 

0.74314483 

-0.06059459 

0.73907379 

0.66546503 

1.02985173 

13 

0.80901697 

0.58778524 

-0.61303794 

0.86128121 

0.49726094 

1.27005431 

14 

0.91354544 

0.40673663 

-1.08526954 

0.94584652 

0.30732416 

1.44404624 

15 

0.97814760 

0.20791169 

-1.35587338 

0.98907379 

0.10395584 

1.53488547 

16 

0.99999999 

0. 

-1.35587353 

0.98907379 

-0.10395584 

1.53488551 

17 

0.97814760 

-0.20791169 

-1.00527014 

0.94584651 

-0. 30732416 

1.44404645 

13 

0.91354544 

-0.40673663 

-0.61303915 

0.86128121 

-0.49726094 

1.27005479 

19 

0.80901699 

-0.58778524 

* 

-0.06059557 

0. 73907379 

-0.66546503 

1.02985221 

20 

0.66913059 

-0.74314483 

0.43530708 

0. 58456530 

-0.8045H51 1 

0.75146053 

21 

0.50000000 

-0.86602540 

0.77778012 

0.40450849 

-0.90854095 

0.47140203 

22 

0.30901697 

-0.95105650 

0.95146379 

0.20677272 

-0.97278919 

0.22030934 

23 

0. 10452846 

-0.99452189 

-0. 

-0.99452189 

0.00000038 

0.99999999 

24 

-0. 10452846 

-0.99452189 

-0.20677273 

-0.972*78919 

-0.22030871 

0.95146407 

25 

-0. 30901700 

-0.95105650 

-0.40450849 

-0.90854094 

-0.47140150 

0.77778063 

26 

-0.50000000 

-0.86602539 

-0.58456530 

-0.80458509 

-0.75146005 

0.43530780 

27 

-0.6691 3061 

-0. 74314480 

-0. 73907380 

-0.66546501 

-1.02985184 

-0.06059480 

28 

-0.80901700 

-0.58778523 

-0.86128122 

-0.49726093 

-1.27005456 

-0.61303858 

29 

-0.91354544 

-0.40673663 

-0.94584651 

-0. 30732416 

-1.44404650 

-1.08527029 

30 

-0.97814760 

-0.20791168 

-0.98907375 

-0. 10395583 

-1.5  3480554 

-1.35587361 

31 

-0.  19999990 

0.00000001 

SIGMA 

0.18791410 
0.54974247 
0.86980452 
1.12206495 
1.28728521 
1.36378153 
1.37810113 
1.37502626 
1.37810123 
1.36378188 
1.28728572 
1.12206513 
0.86980546 
0.54974321 
0.18791476 
-0.18791421 
-0.54974271 
-0.86980455 
-1.12206489 
-1.28728540 
-1.36378179 
-1.3781011  1 
-1.37502635 
-1.37810141 
-1.36378205 
-1.28728585 
-1.12206538 
-0.86980533 
-0.54974335 
-0.18791355 


Figure  14.  -  Continued 


(c)  Solution  at  90°  angle  of  attack 
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DOUGLAS  AIRCRAFT  COMPANY 
LONG  BEACH  DIVISION 


2-0  CASCAOb  TEST  PROHLtM 
CASE  ONE 

NON-UN I  FORM  UNSFT  FLOW  SOLUTION  NO.  1 
UN TRANS FORM ED  COORDINATES 


X 

Y 

V 

CP 

1 

-0. 79999999 

-0.00000001 

-0.98907378 

0.  10395583 

0.25337435 

0.93580144 

2 

-0.97814759 

0.20791 168 

-0.94584651 

0.30732416 

0.23853764 

0.94309979 

3 

-0.91354544 

0.40673663 

-0.86128122 

0.49726093 

0.21024562 

0.95579678 

4 

-0.e0901700 

0.58778523 

-0. 73907379 

0.66546501 

0.17160053 

0.97055326 

5 

-0.66913059 

0.74314480 

-0.58456529 

0.80458509 

0. 12796326 

0.98362540 

6 

-0.49999999 

0.86602539 

7 

-0.40450649 

-0.30901700 

0.90854094 

0.9^105650 

0.08709732 

0.99241406 

-0.20677274 

0.97278919 

0.05785099 

0.99665326 

8 

-0. 10452846 

0.99452189 

0. 

0.99452189 

0.04723441 

0.99776891 

9 

0.  10452846 

0.99452189 

0.20677272 

0.97278919 

0.05785099 

0.99665326 

10 

0.30901699 

0.95105650 

0.40450849 

0.90854095 

0.08709731 

0.99241406 

11 

0.5o000000 

0.86602540 

0.68456530 

0. 804585*1 1 

0.12796327 

0.98362540 

12 

0.66913061 

0.74314483 

0.73907379 

0. 66546503 

0. 17160054 

0.97055325 

13 

0.80901697 

0.58778524 

0.86128121 

0.49726094 

0.21024559 

0.95579679 

14 

0.91354544 

0.40673663 

0.94584652 

0.30732416 

0.23853765 

0.94309978 

15 

0.97814760 

0.20791 169 

0.98907379 

0.10395584 

0.25337436 

0.93580143 

16 

0.99999999 

0. 

0.98907379 

-0. 10395584 

0.25337438 

0.93580142 

17 

0.97814760 

-0.20791 169 

0.94584651 

-0.30732416 

0.23853766 

0.94309978 

16 

0.91354544 

-0.40673663 

0.86128121 

-0.49726094 

0.21024563 

0.95579677 

19 

0.80901699 

-0.58778524 

0.73907379 

-0.66546503 

0.17160054 

0.97055325 

20 

0.66913059 

-0.74314483 

0.58456530 

-0.80458511 

0. 12796328 

0.98362540 

21 

0.50000000 

-0.86602540 

0.40450849 

-0.90854095 

0.08709731 

0.99241406 

22 

0.30901697 

-0.95105650 

0.20677272 

-0.97278919 

0.05785100 

0.99665326 

23 

0.  10452846 

-0.99452189 

-0. 

-0.99452189 

0.04723441 

0.99776891 

24 

-0.  10452846 

-0.99452 189 

-0.20677273 

-0.97278919 

0.05785100 

0.99665326 

25 

-0.30901700 

-0.95105650 

-0.40450849 

-0.90854094 

0.08709732 

0.99241405 

26 

-0.50000000 

-0.86602539 

-0.58456530 

-0.80458509 

0.12796328 

0.98362540 

27 

-0.66913061 

-0.74314480 

-0. T3907380 

-0.66546501 

0.17160055 

0.97055325 

28 

-0.80901700 

-0.58778523 

-0.86128122 

-0.49726093 

0.21024563 

0.95579677 

29 

-0.91354544 

-0.40673663 

-0.94584651 

-0.30732416 

0.23853769 

0.94309977 

30 

-0.97814760 

-0.20791168 

-0.98907375 

-0.10395583 

0.25337437 

0.93580142 

31 

-C. 99999990 

0.00000001 

SIGHA 

-0.02108791 
-0.06055873 
-0.09188694 
-0.10970391 
-0.10935894 
-0.08856266 
-0.04960032 
-0.00000000 
0.04960030 
0.03856265 
0.10935894 
0.10970389 
0.09188698 
0.06055875 
0.02108794 
-0.02108793 
-0.06055875 
-0.09188694 
-0.10970390 
-0. 10935894 
-0.08856265 
-0.04960031 
-0.00000000 
0.04960032 
0.08856265 
0.10935895 
0. 10970390 
0.09188696 
0.06055877 
0.02108780 


Figure  14.-  Continued 


(d)  Circulatory  flow  solution 
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2-0  CASCADE  TEST  PROBLEM 


SPACING  * 

3.00000000 

ALPHA  « 

9.99999998 

DELTA  ALPHA  »■ 

19.60917735 

INLET  ALPHA  * 

19.51188397 

V  INLET  * 

1.04480879 

XHC  * 

0. 

EXIT  ALPHA  * 

-0.09729340 

V  EXIT  = 

0.98480917 

YMC  « 

0. 

COMBINED  VELOCITIES 

BODY  NU.  1  UNTRANSFGRHED 

COORDINATES 

X 

Y  - 

VC 

CP 

DELTA 

S 

1 

-0.99999999 

-0.98907378 

-O.OOOOOCOl 

0.10395583 

-0.27208283 

0.92597093 

0.20905692 

2 

-0.97814759 

-0.94584651 

0.20791168 

0.30732416 

-0.81635635 

0.33356231 

0.20905693 

3 

-0.91354544 

-0.86128122 

0.40673663 

0.49726093 

-1.36190663 

-0. 854  7d967 

0.20905690 

4 

-0.80901700 

-0.73907379 

0.58778623 

0.66546501 

-1.90977612 

-2.64724481 

0.20905694 

5 

-0.66913059 

-0.58456529 

0.74314480 

0.80458509 

-2.4524780? 

-5.01464844 

0.20905693 

6 

-0.49999999 

-0.40450849 

0.86602539 

0.90854094 

-2.95312697 

-7.72095883 

0.20905690 

7 

-0.30901700 

-0.20677274 

0.95105650 

0.97278919 

-3.32314461 

-10.04329002 

0.20905693 

8 

-0.10452846 

0. 

0.99452189 
0.99452  189 

-3.44250140 

-10.85081577 

0.20905694 

9 

0.10452846 

0.20677272 

0.99452189 

0.97278919 

-3.24663243 

-9.54062212 

0.20905691 

10 

0.30901699 

0.40450849 

0.95105650 

0.90854095 

-2.78941125 

-6.78081506 

0.20905693 

11 

0.50000000 

0.58456530 

0.86602540 

0.80458511 

-2. 19149938 

-3.80266953 

0.20905692 

12 

0.66913061 

0.73907379 

0.74314483 

0.66546503 

-1.55211282 

-1.40905419 

0.20905691 

13 

0.80901697 

0.86128121 

0.58778524 

0.49726094 

-0.92082220 

0.15208647 

0.20905692 

14 

0.91354544 

0.94584652 

0.40673663 

0.30732416 

-0.31484519 

0.90087250 

0.20905692 

15 

0.97314760 

0.98907379 

0.20791 169 
0.10395584 

0.26097684 

0.93189109 

0.20905692 

16 

0.99999999 

0.98907379 

0. 

-0. 10395584 

0.80514307 

0.351 74464 

0.20905692 

17 

0.97814760 

0.94584651 

-0.20791 169 
-0. 30732416 

1.31820282 

-0.73765868 

0.20905692 

18 

C. 91354544 
0.86128121 

-0.40673663 

-0.49726094 

1.80423087 

-4 .25524902 

0.20905691 

19 

0.80901699 
0. 73907379 

-0.58778524 

-0.66546503 

2.27079713 

-4. I56r 1959 

0.20905694 

20 

0.6691 3059 
0.58456530 

-0.74314483 

-0.80458511 

2.72169322 

-6.40761393 

0.20905691 

21 

0.50000000 

0.40450849 

-0.H6602540 

-0.90854095 

3.13636646 

-8.83679450 

0.20905694 

22 

0.30901697 

0.20677272 

-0.95105650 

-0.97278919 

1.44485435 

-10.86/02144 

0.20905690 

23 

0.  10452846 
-0. 

-0.99452189 

-0.99452189 

3.54187548 

-11.54488182 

0.20905694 

24 

-0. 10452846 
-0.20677273 

-0.99452189 

-0.97278919 

3.36834201 

-10.34572780 

0.20905692 

25 

T0. 30901700 
-0.40450849 

-0.95105650 

-0.90854094 

2.97265065 

-7.83665180 

0.20905691 

26 

-0.50000000 

-0.58456530 

-0.86602539 

-0.80458509 

2.46071422 

-5.05511445 

0.20905693 

27 

-0.6691 3061 
-0.73907380 

-0.74314480 

-0.66546501 

1.91313393 

-2.66008145 

0.20905692 

28 

-0.80901700 
-C. 86128122 
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Figure  14.-  Continued 
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Figure  14.  -  Continued  (f )  Combined  solutions  for  off-body  points 
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Figure  15.-  Continued 
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Att :  Librarian 

North  Carolina  State  College 
Head  of  liechanical  Engineering  Dept. 
Raleigh,  North  Carolina 
Att:  Dr.  R.  W.  Truitt 

Ohio  State  University 
Research  Foundation 
Columbus  10,  Ohio 
Att:  Library 


Pennsylvania  State  University 
University  Park,  Pennsylvania 
Att:  Professor  George  F.  Wislicenus 

Purdue  University 

School  of  Aeronautical  Engineering 

Lafayette,  Indiana 

Att:  Library 

University  of  Southern  California 

Engineering  Center 

935  West  37th  Street 

Los  Angeles  Y}  California 

Att:  Dr.  Raymond  Chuan 

Stanford  University 
Department  of  Me  thematic  s 
Stanford,  California 
Att:  Dr.  B.  Perry 

Stanford  University 
Department  of  Mathematics 
Stanford,  California 
Att:  Dr.  E.  Y.  Hsu 

Stanford  University 
Department  of  Mathematics 
Stanford,  California 
Att:  Profe  I.  Flugge-Lotz 

Stevens  Institute  of  Technology 
5th  and  Hudson  Street 
Hoboken,  New  Jersey 
Att:  Dr.  John  P.  Breslin 

Stevens  Institute  of  Technology 

Hoboken,  New  Jersey 

Att:  Director,  Davidson  Laboratory 

Stevens  Institute  of  Technology 
Davidson  Laboratory 
Hoboken,  New  Jersey 
Att:  Mr.  C.  J.  Henry 

Stevens  Institute  of  Technology 
Davidson  Laboratory 
Hoboken,  New  Jersey 
Att:  Mr.  S.  Tsakonas 

University  of  Texas 
Defense  Research  Laboratory 
P.  0.  Box  8029 
Austin,  Texas 

Att:  Professor  M.  J.  Thompson 


University  of  Washington 

Dept  of  Aeronautical  Engineering 

Seattle  5,  Washington 

Att:  Professor  R.  E.  Street 

State  University  of  Iowa 

Iowa  Institute  of  Hydraulic  Research 

Iowa  City,  Iowa 

Att:  Prof.  L.  Landweber 

State  University  of  Iowa 

Iowa  Institute  of  Hydraulic  Research 

Iowa  City,  Iowa 

Att:  Prof.  Hunter  Rouse 

University  of  California 
Dept  of  Naval  Architecture 
Berkeley,  California 
Att:  Prof.  H.  A.  Sc’hade,  Head 

St.  Anthony  F  11s  Hydraulic  Laboratory 
University  of  Minnesota 
Minneapolic,  Minn. 

Att:  Prof,  B.  Silberman 

St.  Anthony  Falls  Hydraulic  Laboratory 
University  of  Minnesota 
Minneapolis,  Minn. 

Att:  Mr.  J.  N.  Wetzel 

Ordnance  Research  Laboratory 
Pennsylvania  State  University 
University  Park,  Penn. 

Att:  Dr.  M.  Sevik 

Virginia  Polytechnic  Institute 
Aerospace  Engineering  Dept. 

Blacksburg,  Virginia 
Att:  Dr.  H.  A.  Hassan 


INDUSTRIAL  AND  RESEARCH  COMPANIES 

Midwest  Research  Institite 
425  Volker  Blvd. 

Kansas  City  10,  Missouri 
Att :  Mr .  Zeydel 

Director,  Dept,  of  Mechanical  Sciences 
Southwest  Research  Institute 
85OO  Culebra  Road 
San  Antonio  6,  Texas 


James  Forrestal  Research  Center 
Dept,  of  Aeronautical  Engineering 
Princeton,  New  Jersey 
Att:  Professor  C.  D.  Perkins,  Chairman 

Department  of  Mechanical  Sciences 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio  6,  Texas 
Att:  Dr.  H.  N.  Abramson 

Department  of  Mechanical  Sciences 
Southwest  Research  Institute 
85OO  Culebra  Road 
San  Antonio  6,  Texas 
Att:  Mr.  G.  Ransleben 

Department  of  Mechanical  Sciences 
Southwest  Research  Institute 
85OO  Culebra  Road 
San  Antonio  6,  Texas 
Att:  Editor,  Applied  Mechanics 
Review 

Boeing  Airplane  Co.,  Seattle  Division 
P.  0.  Box  3107 
Seattle  l4,  Washington 
Att:  Mr.  G.  Schairer 

Boeing  Airplane  Co.,  Seattle  Division 
P.  0.  Box  3107 
Seattle  l4,  Washington 
Att:  Mr.  M.  J.  Turner 

Convair 
P.  0.  Box  1950 
San  Diego  12,  California 
Att:  Mr.  A.  D.  MacLellan,  Systems 
Dynamic  Group 

Convair 
P.  0.  Box  1950 
San  Diego  12,  California 
Att:  Mr.  H.  T.  Brooke,  .Hydrodynamics 
Group 

Electric  Boat  Division 
General  Dynamics  Corp. 

Groton,  Conn. 

Att:  Mr.  Robert  McCandliss 


Gibbs  and  Cox,  Inc. 

21  West  Street 
New  York,  New  York 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,  New  York 
Att:  Mr.  E.  Baird 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,  New  York 
Att:  Mr.  E.  Bower 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,  New  York 
Att:  Mr.  14.  C.  Tilgner 

Grumman  Aircraft  Engineering  Corp. 
Dynamic  Developments  Division 
Babylon,  New  York 

President,  Hydronautics,  Inc. 
Pindell  School  Road 
Laurel,  Howard  County, Maryland 
(2  copies) 

Lockheed  Aircraft  Corp. 

Missiles  and  Space  Division 
Palo  Alto,  California 
Att:  R.  W.  Kermeen 

Lockheed  Aircraft  Corp. 

Missiles  and  Space  Division 
Palo  Alto,  California 
Att:  L.  A.  Rodert 

Technical  Research  Group,  Inc. 
Route  110 

Melville,  Long  Island,  New  York 
Att?  Dr.  Jack  Kotik  (2  copies) 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Att:  Dr.  B.  Parkin 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Atfc:  D.  Morriss 


General  Applied  Sciences  Laboratories,  Inc. 

Merrick  and  Stewart  Avenues  Armour  Research  Foundation 

Westbury,  Long  Island,  New  York  West  35th  Street 

Att:  Dr.  F.  Lane  Chicago  lo,  Illinois 

Att :  Document  Librarian 


Avc o-Eve rett  Research  Laboratory 
2385  Revere  Beach  Parkway 
Everett  49,  Massachusetts 
Att:  Dr.  Richard  H.  Levy 


Lockheed  Missile  Systems  Division 
Research  and  Development  Laboratory 
Sunnyvale,  California 
Att:  Dr.  W.  Griffith 


Avco  Manufacturing  Corp. 
2385  Revere  Beach  Parkway 
Everett  49,  Massachusetts 
Att:  Dr.  A.  Kantrowitz 


Eorair 

Division  of  Northrop  Aircraft,  Inc. 
Hawthorne,  California 
Att:  Dr.  W.  Pfenninger 


AVCO  Research  Library 
2385  Revere  Beach  Parkway 
Everett  49,  Massachusetts 


Mr.  K.  E.  Van  Every 

Norair,  Division  of  Northrop  Aircraft,  Inc. 
Hawthorne,  California 


Aeronutronics  Systems,  Inc 
1234  Air  Way 
Glendale,  California 
Att:  Dr.  L.  Kavanau 

Chance-Vought  Corporation 

Dallas,  Texas 

Att:  W.  C.  Schoolfield 


Sulzer  Bros.,  Inc. 

19  Rector  Street 
New  York  6,  New  York 

North  American  Aviation,  Inc 
4300  East  5th  Avenue 
Columbus,  Ohio 
Att:  Mr.  R.  M.  Crone 


General  Dynamics  Corporation 
Electric  Boat  Division 
Groton,  Connecticut 
Att:  Frank  Walker  Caldwell 

General  Electric  Company 
2900  Campbell  Avenue 
Schenectady  6,  New  York 
Att :  Library 

General  Electric  Company 
Research  Laboratory 
Schenectady,  New  York 
Att:  Dr.  H.  T.  Nagamatsu 

Grumman  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,  New  York 
Att:  Fluid  Mechanics  Section 

Hughes  Aircraft  Company 
Florence  Ave.  ac  Teal  Street 
Culver  City,  California 
Att:  Dr.  A.  E.  Puckett 

United  Aircraft  Corporation 
East  Hartford,  Connecticut 
Att:  Mr.  J.  G.  Lee 

McDonnell  Aircraft  Corporation 
St.  Louis,  Missouri 
Att:  Library 


North  American  Aviation,  Inc. 
Aerophysics  Department 
12214  Lakewood  Boulevard 
Downey,  California 
A.tt:  Dr.  E.  R.  van  Driest 

Republic  Aviation  Corporation 
Conklin  Street 

Farmingdale,  Long  Island,  New  York 
Att:  Dr.  W.  J.  O'Donnell 

Therm  Incorporated 
Aerophysics  Section 
Ithaca,  New  York 
Att:  D.  E.  Ordway,  Head 


MISCELLANEOUS 

Pacific  Aeronautical  Library 
Institute  of  the  Aeronautical 
Sciences 

7660  Beverly  Boulevard 
Los  Angeles  36,  California 

Applied  Mechanics  Review 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas 


u 


Society  of  Naval  Architects  and 
Marine  Engineers 
7^  Trinity  Place 
New  York,  New  York 

Yon  K.  Jacob 
Computing  Center 

California  Institute  of  Technology 
1201  East  California  Street 
Pasadena  3,  California 


u 


